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Abstract. In this paper we continue investigation of the constant astigmatism equation 
z yy + (l/z) xx + 2 = 0. Wc newly interpret its solutions as describing spherical orthogonal 
equiareal patterns, with relevance to two-dimensional plasticity. We show how the classical 
Bianchi superposition principle for the sine-Gordon equation can be extended to generate 
an arbitrary number of solutions of the constant astigmatism equation by algebraic 
manipulations. As a by-product, we show that sine-Gordon solutions give slip line fields on 
the sphere. Finally, we compute the solutions corresponding to classical Lipschitz surfaces 
of constant astigmatism via the corresponding equiareal patterns. 



1. Introduction 

It is well known that the classical Backlund transformation [4] for the sine-Gordon 
equation u^ v = sin-u as well as the Bianchi permutability property [7] have been 
discovered in the context of pseudospherical surfaces, i.e., surfaces of constant negative 
Gaussian curvature. It is perhaps less known that historical roots of these developments 
lie in another class of surfaces, characterised by the constancy of the difference p 2 — pi 
between the principal radii of curvature pi, P2', see [26] for the historical account. 
Lying covered with dust and oblivion for almost a century, the surfaces satisfying 
P2 — Pi = const reemerged recently from the systematic search for integrable classes 
of Weingarten surfaces conducted by Baran and one of us [5] . Although nameless in the 
nineteenth century, in [5] they have been named the surfaces of constant astigmatism 
in connotation with the astigmatic interval [31] of the geometric optics, albeit without 
suggesting any specific application. 

Undoubtedly, the most important results about constant astigmatism surfaces are 
due to Bianchi. In [6] (see also [8, §130]), Bianchi observed that evolutes (i.e., focal 
surfaces) of surfaces satisfying p-i — p\ = const are pseudospherical. In the same paper 
he also constructed surfaces satisfying p 2 — p\ = const as involutes corresponding to 
parabolic geodesic systems on pseudospherical surfaces. Apparently, Bianchi was the 
first to obtain surfaces of constant astigmatism explicitly, namely, surfaces [6, eq. (30)] 
corresponding to Dini's pseudospherical helicoids (see, e.g., [27, §1.4.2] or [32, p. 183]). 
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Lipschitz [22] obtained another class of surfaces of constant astigmatism; within the 
full class given in terms of elliptic integrals he pointed out a subclass of surfaces of 
revolution, further investigated by von Lilienthal. 

Let us stress that the aforementioned constructions of Lipschitz and Bianchi refer 
to ad hoc parameterisations. Bianchi used the rotation angle and a parameterisation of 
the generating tractrix of the helicoid, while Lipschitz employed spherical coordinates 
on the Gaussian sphere. In [5] we observed that under an adapted parameterisation 
by lines of curvature the constant astigmatism surfaces correspond to solutions of the 
constant astigmatism equation 



(x,y are natural parameters in the sense of Ganchev and Mihova [16]). The geometric 
link to pseudospherical surfaces induces a nonlocal transformation to the sine-Gordon 
equation and vice versa. Since curvature coordinates on constant astigmatism surfaces 
correspond to parabolic geodesic coordinates on the pseudospherical surfaces, and these 
are not the coordinates the sine-Gordon equation is referred to, the transformations 
change both the dependent and independent variables. Explicit formulas can be found 
in [5], ready to be applied to the sine-Gordon solutions, which are known in abundance, 
see [1, 14, 24, 25] and references therein. However, the only explicit instance of such a 
relationship we were able to find in the literature was that of the Bianchi surfaces [6, 
eq. (30)] to the Dini helicoid; Fig. 1 presents a plot of them as unparameterised surfaces. 



In this paper we continue the investigation of the constant astigmatism equation (1) 
and its solutions. Firstly, we show that equation (1) describes orthogonal equiareal 
patterns (Sadowsky [29, 30]) on the sphere, i.e., a system of local coordinates x*,x 2 
such that the metric coefficients satisfy g^ = 0, det g = 1. Hence, the area element is 
simply da; 1 A dx 2 and the area of the curvilinear rectangle a 1 < x l < b l , % — 1, 2, is equal 
to (b 1 — a l ){b 2 — a 2 ). It follows that the curvilinear rectangles formed by "uniformly 
spaced" coordinate lines are of equal area, which explains the terminology. 




(1) 




Figure 1. Dini's pseudospherical surface (left) and its constant astigmatism involute (right) 
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The equiareal property in plasticity theory can be traced back to Boussinesq [11]. 
Seventy years later, Sadowsky [29, 30], rediscovered the "equiareal patterns" as 
configurations of the principal stress lines under the Tresca yield condition (see [19]) 
and gave them their name. Later Hill [18] gave a kinematic interpretation of these 
patterns. Coburn [13, Thm. 1] established the same equiareal property, this time for 
slip lines under a different yield condition. Ament [2] discovered a relation to the 
class of Weingarten surfaces, determined by the constancy of the difference between 
the principal curvatures (as opposed to the difference between the principal radii of 
curvature). Finally, Fialkow [15, Th. 4.1] observed relevance of orthogonal equiareal 
patterns to conformal geometry. 

The contents of this paper are as follows. Section 2 contains the necessary 
background. In Section 3, we observe that adapted curvature coordinates on constant 
astigmatism surfaces correspond to orthogonal equiareal patterns on the Gaussian 
sphere. Inspired by the aforementioned relation to plasticity, we construct a two- 
dimensional stress tensor formally satisfying both the Tresca yield condition and the 
equilibrium equations. Physical relevance of our purely mathematical construction is 
not a primary concern, yet the flow of a thin plastic layer around a sphere seems 
to be a realistic picture. Guided by this picture we investigate the maximum shear 
stress directions, positioned at the angle of 7r/4 to the principal stress directions. The 
corresponding trajectories are known as slip lines; we show them to be related to 
solutions of the sine-Gordon equation later in section 4. 

Section 4 is devoted to a simplified reconstruction of constant astigmatism surfaces 
from a pair of complementary pseudospherical surfaces [6] or [8, §136], i.e., under 
the frequently occurring condition that both evolutes are known. Complementary 
pseudospherical surfaces are easy to find among those resulting from the famous 
and powerful Bianchi permutability theorem [7]. It turns out that given a pair of 
complementary pseudospherical surfaces, the corresponding (unparameterised) constant 
astigmatism surface can be obtained by pure algebraic manipulations and differentiation. 
This is also true for geodesies (as proved by Bianchi himself) and, hence, for one of 
the curvature coordinates, while obtaining the other requires one integration. However, 
owing to a suitable extension of the Bianchi superposition principle this integration 
needs to be done only once. We also observe (Proposition 4) that the coordinates £, i] 
the sine-Gordon equation is referred to correspond to slip line fields on the spherical 
image of the constant astigmatism surface. 

Finally, in section 5 we pay another longstanding debt and find the function z(x, y) 
corresponding to the Lipschitz surfaces. As we already mentioned, Lipschitz computed a 
class of constant astigmatism surfaces in terms of spherical coordinates on the Gaussian 
image. The result being not easily transformable to curvature coordinates, we compute 
the associated orthogonal equiareal pattern directly from the definitions to observe that 
solutions of the Lipschitz class are invariant solutions with respect to Lie symmetries. 
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2. Preliminaries 

In this section we recall previous results about the constant astigmatism surfaces; see [5] 
for details. We consider surfaces immersed in Euclidean space under parameterisation 
by the lines of curvature (also known as curvature coordinates). Hence, the fundamental 
forms can be written as 

I = u 2 dx 2 + v 2 dy 2 , II = - dx 2 + - dy 2 , III = t> dx 2 + ^dy 2 , 

Pi P2 Pi P 2 

where pi, p% are the principal radii of curvature. The first two forms determine the 
surface up to the rigid motions (Bonnet theorem). 

A surface is said to be of constant astigmatism if the difference p2 — pi between the 
principal radii of curvature is a nonzero constant (if zero, then the surface is a part of 
the sphere). We assume the ambient space to be scaled so that p 2 — Pi — ±1- 

Definition 1. A parameterisation by lines of curvature is said to be adapted if 

uv(—-—)=±l (2) 

V Pl p2 J 

holds. 

This is the natural parameterisation recently introduced by Ganchev and 
Mihova [16, Prop. 5.6] with the arbitrary constant being normalised to ±1. Every 
constant astigmatism (more generally, Weingarten) surface can be equipped with an 
adapted parameterisation by lines of curvature, see [16, Prop. 5.6] or [5]. Henceforth we 
assume that x, y are adapted coordinates. Then, according to [5], the nonzero coefficients 
of the three fundamental forms of a surface of constant astigmatism can be expressed 
through a single variable z(x,y): 

z 5 (lnz — 2) In ,2 lnz — 2 Inz 

u = , v = ^' P l = o ' P2 = ^. 

Z 2z 2 * * 

Obviously, condition (2) is satisfied. 

Let r(x,y) be the surface of constant astigmatism corresponding to z(x,y), let 
n(x, y) denote the unit normal vector. Then r, n satisfy the Gauss- Weingarten system 

(In z)z x (\nz — 2)zz v 1 „ x 

(In z — 2)zz x 

(3) 



2{\nz- 


-2)/ : 


(lnz) 


Zy 

r. 


2(mz - 


-2)z ■ 


(lnz 


)z x 



Tx v-<^-~ ^~ Tx 2\nz ryi 



(lnz — 2)z y lnz 



Vyy ~ 2{\nz-2)z* Vx " 2z\nz Yy+ 2z n ' 

2 2 
nx ~~\nz^2 rx ' ny ~~\nz Vy - 
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Note that ei = r x /u, e 2 = r y/v, and n = ei x e 2 constitute an orthonormal frame. 

Compatibility conditions of the Gauss-Weingarten system constitute the Gauss- 
Mainardi-Codazzi system, which in our case amounts to the Gauss equation alone, and 
coincides with the constant astigmatism equation (1). 

According to Bianchi [6] (see also [9, §136]), if r is a surface of constant astigmatism 
and n is its normal, then the two evolutes 

r + pin, r + p 2 n (4) 

are pseudospherical surfaces. These are said to be complementary. 

For further reference, we also recall a list of symmetries of equation (1). Lie 
symmetries are completely known, see [5]. They are the x-translation %^(x,y,z) = 
(x + c,y,z), the y-translation % y c {x,y,z) = (x,y + c,z), and the scaling & c (x,y,z) = 
(e~ c x,e c y,e 2c z), where c is a real parameter. 

We shall also refer to a discrete symmetry 3(x,y,z) = (y,x,l/z), called the 
involution. Obviously, 

3 o 3 = Id, 

3o%z = %y a o3, 

& c o%: = %: /c o 6 C , 

& c o 3 = 3 o 6i /c . 

Translations are mere reparameterisations of the corresponding constant astigmatism 
surface. The scaling symmetry corresponds to an offsetting, i.e., takes a surface to a 
parallel surface (moves every point a unit distance along the normal). The involution 
interchanges x and y (swaps the orientation), followed by a unit offsetting. 



3o%Z = ZZo3, 
6 C o S» = o 6 C , 



3. Orthogonal equiareal patterns and slip line fields 

The geometric meaning of the variable z can be seen from the third fundamental form, 
which turns out to be simply 

III = zdx 2 + l dy 2 . 

z 

Since III = dn ■ dn coincides with the first fundamental form of the Gaussian sphere 
n(x, y), it follows that one obtains a rather special parameterisation of the latter. 

Definition 2. By an orthogonal equiareal pattern on a surface S we shall mean a 
parameterization x, y such that the corresponding first fundamental form is 

l s = zdx 2 + 1 dy 2 , (5) 

z 

z being an arbitrary function of x, y. 
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Let R denote the position vector of a point on the surface S. Since detls = 1, the 
local parameterisation H(x, y) is an area preserving map from the plane to the surface 
S. Moreover, the coordinate lines are, obviously, orthogonal. These two properties imply 
that evenly distributed coordinate lines cover the surface with curvilinear rectangles of 
equal area (see the Introduction). 

Example 1. The Archimedean projection. A simple example of an orthogonal equiareal 
pattern on the sphere that can be seen on Fig. 2 is delivered by the well-known 
Archimedean projection of the cylinder (cos?/, sin y, x) onto an inscribed sphere. In this 
case, (x,y) is sent to (y/l — x 2 cosy, Vl — x 2 siny, x) and we have 

dr 2 

iArch = ~ 2 )dy , 

1 — X 

i.e., z = 1/(1 — x 2 ). According to [5], this solution of the constant astigmatism equation 
corresponds to von Lilienthal surfaces [21]. 



Figure 2. The Archimedean equiareal parameterisation of the sphere 



Not only every constant astigmatism surface generates an orthogonal equiareal 
parameterization of the unit sphere; a converse statement is also available. 

Proposition 1. Let n(x,y), ||n|| = I, be an orthogonal equiareal pattern on the unit 
sphere S. Then z defined by formula (5) is a solution of the constant astigmatism 
equation (1). 

Proof. Using the well-known Brioschi formula to compute the Gaussian curvature of the 
sphere, we obtain 

2 Zyy 2\z) xx 

The constant astigmatism equation (1) easily follows. □ 
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The corresponding constant astigmatism surface can be reconstructed from the last 
two equations of the Gauss- Weingarten system (3). 

Let us stress that all the point symmetries given in Sect. 2 can be understood as 
reparameterisations of the corresponding orthogonal equiareal pattern on the Gaussian 
sphere. In particular, scaling & c means shrinking the pattern along one family of lines, 
compensated by stretching it along the orthogonal family of lines. 

In the case of S being a plane, the notion of an orthogonal equiareal pattern was 
introduced by Sadowski [29, 30] in the context of two-dimensional plasticity. Choosing 
the vectors d x ,d y along the principal stress directions (i.e., eigenvectors of the stress 
tensor op, Sadowski derived the equiareal property from the equilibrium condition 
div a — and the Tresca yield condition o\ — a\ — const. 

Let us reverse the line of reasoning and reconstruct a two-dimensional stress tensor 
from a given orthogonal equiareal pattern g = Is- In what follows, all components are 
taken with respect to the basis d x ,d y of the tangent space and indices are raised and 
lowered with the metric. 

Proposition 2. Consider an orthogonal equiareal pattern g = gij dx % dx J such that 

9n = z, gu = 92i = 0, g 22 = l/z. 
Then the tensor a given by the components 

a\ = \ In 2, o\ = o\ = 0, al = \{\nz- 2). (6) 
satisfies <7?j = (the equilibrium equation) and o\ — o\ = \ (the Tresca yield condition). 
Proof. From the metric coefficients we produce the Christoffel symbols 

pi _ _p2 _ %x_ pi _ _p2 _ Zy_ pi _ Z x p2 _ zz y 
111- i 12- 2 ^' L 12~ L 22~ 2z , L 22~ 2z 3> 1 11 ~ 3 

as well as the twice contravariant tensor 

<7 U = — , a 12 = a 21 = 0, a 22 = h ^ 2 z. (7) 
2z' '2 K ' 

The yield condition o\ — o\ = 1 is obvious, while checking the equilibrium equation 
a 1 ? = is a matter of routine. □ 

This proposition holds for any surface S equipped with an orthogonal equiareal 
pattern g. Note that if S is a unit sphere, then the Tresca yield condition follows from 
the constant astigmatism property, since a\ = p 2 and a\ = p\. Conversely, if o~\ and 
o\ are arbitrary functions of z, then the equilibrium equation and the yield condition 
imply the same a\ and a\ as in (6) up to an additive and a multiplicative constant. 

In the rest of this section we recall the derivation of the Mohr circle and slip lines 
in two-dimensional plasticity (see, e.g., [19]). We consider a symmetric stress tensor a 
diagonalised along the principal stress directions d x ,d y , i.e., 

el =P, o\ = o\ = 0, a\ = q, g 12 = 0. 
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Let w be a unit vector with components (cos 0/v / gn , sin0/\/^22 ) with respect to the 
basis d x ,d y . The stress er(w) = pcos0/v / gn + 1 sin / 'V922 can be decomposed into a 
sum of the normal stress ctjv(w) and the shear stress o"t(w), where <tat(w) _L <7t(w) 
and, by definition, oat(w) = aw is a multiple of w. Hence, 

cr T (w) = cr(w) - aw = ((p - a) cos^/V^n , (g - a) sin 0/ ) . 
Obviously, ||<jjv(w)|| = |a|; likewise, we introduce f3 = ||<7 T (w)||. By orthogonality, 

= <7jv(w) • (j T (w) = a(p — a) cos 2 + a(g — a) sin 2 0. 
Excluding cos and sin from the last equation and the condition 

= (3 2 — ot(w) • ctt(w) = (p — a) 2 cos 2 + (g — a) 2 sin 2 0, 
we conclude that all admissible values of a, j3 belong to the Mohr circle [23] 

It follows that the extremal values a = p, q of the normal stress magnitude a are 
achieved when /3 — 0, i.e., when w lies in one of the principal stress directions, as it 
should be. The extremal values \\{p — q)\ of the shear stress magnitude (5 are achieved 
when a — |(p + q). To find the corresponding vectors w, we determine the acute angle 
between w and ±d x = (±1, 0). Substituting a = \{p + q) into cos 2 = {ot — q)/(p — q), 
we obtain cos 2 = |, meaning that = \ix. 

Now, the Tresca criterion (see, e.g., [19]) says that yielding occurs whenever the 
maximal shear stress magnitude (3 achieves a threshold depending on the material. It 
follows that the stress tensor satisfies p — q = const, which is called the Tresca yield 
condition. The lines along the maximal shear stress direction are called slip lines and, as 
we have already seen, have a constant deviation of \n from the principal stress directions. 

Definition 3. By a slip line field associated with the orthogonal equiareal pattern (5) 
on a surface S we shall mean a parameterization £,1] such that the angle between d x 
and as well as the angle between d y and d v is equal to \-n. 

It follows that slip lines form an orthogonal net. Note the available freedom of 
reparameterisation of each £ or r\ separately. 

Remark 1. In plane plasticity slip line lines form what is called a Hencky net [17]. 
These have been fully described by Caratheodory and Schmidt [12]; a full description 
of orthogonal equiareal patterns in the plane follows. However, the famous Hencky 
conditions satisfied by slip lines in the plane fail on surfaces of non-vanishing curvature. 

Example 2. Continuing Example 1, we easily see that the corresponding orthogonal net 
of slip lines is, by definition, formed by the ±45° loxodromes (lines of constant bearing); 
see Fig. 3 or model No. 249 in the Gottingen collection [28]. Also compare with Zelin's 
superplastic sheet stretched with a spherical punch [33, Fig. 5b]. 

In the next section (Prop. 4) we shall see that solutions of the sine-Gordon equation 
produce slip line fields on the Gaussian sphere of the associated constant astigmatism 
surface. 
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4. Backhand transformation and the superposition principle 

Available already at the beginning of the nineteenth century, the classical Backlund 
transformation [4] in combination with the Bianchi permutability theorem [7] is a 
powerful way to generate pseudospherical surfaces and solutions of the sine-Gordon 
equation. In this section we extend these methods to generate constant astigmatism 
surfaces, solutions of the constant astigmatism equation as well as equiareal patterns 
and slip line fields on the sphere. 

To start with we briefly recall the Backlund transformation [4], see, e.g., [9, 27, 32]. 
Let us consider a pseudospherical surface r(£, 77), where the parameters £,77 are both 
Chebyshev and asymptotic (which is always possible), i.e., 

I = d£ 2 + 2 cos(2w) d£ dr] + dr] 2 , II = 2 sin(2w) d£ drj. 

The position vector r(£, rj) and the unit normal n(£,r?) satisfy the Gauss- Weingarten 
system 

sin 4u 2 
sin 2cu sm2w 
r^ = (sin2o;)n, 

sin Au 2 

^ = sirT^^ " shT2^' (8) 
sin 4u 1 
2 sin 2tu sin 2w 
sin 4u 1 
' 2sin 2 2w " sin2w 6 

The integrability conditions of the above system reduce to the sine-Gordon equation 
uj£ V = I sin 2u. 
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The Backlund transform of our surface is 

+ 1 + A 2 V sin(2a;) ?+ sin(2u;) V' 1 J 

where u/ A ) is another sine-Gordon solution, obtained from the pair of compatible first- 
order equations 



CO 



W = Wf + A sin(w (A) + w), wJ A) = -w„ + ^ sin(w (A) - u). (10) 

Here A is a constant called the Backlund parameter. 

The Backlund transformation is particularly useful in combination with Bianchi's 
permutability theorem [7]; see also, e.g., [9, 27]. To simplify exposition, we shall write 
B^oo to denote a solution u/ A ) of system (10) for a specified value of the integration 
constant c. The Bianchi permutability theorem says that given a pair of Backlund 
parameters Ai ^ A2, then for every choice of integration constants c\,C2 there is a 
unique choice of integration constants c[, c' 2 such that 

B^ 2) Bi Xl) u = B^B^u (11) 

and, moreover, denoting by uj( XiX ^ the common value in (11), then uj( XiX ^ can be 
obtained from the superposition principle 

w (AiA 2 )_ w Ai + A 2 u iXl) -J X2) . . 

tan = — tan (12) 

2 Ai — A2 2 

without further integration, by purely algebraic manipulations. The corresponding 
pseudospherical surfaces r^ AlA2 ^ can be easily obtained by iterating formula (9). 

Assume now that a general solution of system (10) is known for every value of 
the Backlund parameter A. Substituting u = into formula (12), one can also 

compute the Backlund transform w (AlA2As) = £ (A3) w (AlA2) = i3 (A2) w (AlA3) , by purely 
algebraic manipulations. In principle, this process may be repeated indefinitely, leading 
to solutions k;( A i A 2'" A s) depending on any finite number of Backlund parameters and 
integration constants, which are thereby obtained by purely algebraic manipulations. 
Needless to say, the corresponding pseudospherical surfaces r( AlA2 "' As ) can be computed 
by iterating the formula (9). Having summarized the Backlund transformation and 
the Bianchi superposition principle, we proceed to generation of surfaces of constant 
astigmatism by purely algebraic means. 

To start with, we remind the reader that in the particular case of A = 1 
the Backlund transformation coincides with Bianchi's [6] complementarity relation, 
cf. (4) (actually, the Backlund transformation is a combination of complementarity 
and Lie's transformation, and the latter is identity if A = 1). Otherwise said, 
the complementary pseudospherical surfaces result from the particular case B^ of 
the Backlund transformation. Consequently, the superposition formula (12) yields a 
method to obtain abundant pairs of complementary sine-Gordon solutions a/ AlA2 "' As ) 



1 
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and w ( a i a 2-- a s i)^ Likewise, one can also obtain abundant pairs of complementary 
pseudospherical surfaces r^ AlA2 '" A ^ and r ( A i A 2-- A s i) us i n g formula (9). 

Let r be a pseudospherical surface, corresponding to a sine-Gordon solution ou. 
Substituting A = 1 into formulas (9) and (10), we immediately see that the 
complementary surface is 



r 



(i) 



r+ — —. — ^r £ + — — ^r„, (13) 
sin(2w) 5 sin(2w) v ; 



where co^ is the complementary solution of the sine-Gordon equation, satisfying the 
compatible first-order equations 



,W =u i + sin(J 1) = -u;,, + sin(u; (1) -ou). (14) 



Before proceeding further, we recall two important observations due to Bianchi [6] (see 
also [9, §386]). Considering the dependence of cu^ on the integration constant c and 
denoting / = h^du/^/dc), differentiation of (14) gives 

ft = cos(w (1) + to), f v = cos(w (1) - u). (15) 

Similarly, taking one more derivative /' = df/dc, we get 

4 = -e / sin(w (1) +oj), f' v = -e / sin(w (1) - u). (16) 

It follows that, knowing solutions of system (14), we can also obtain solutions of 
systems (15) and (16) by purely algebraic manipulations and differentiation. 

All this is important because surfaces of constant astigmatism are easier to obtain 
from a pair of complementary pseudospherical surfaces r and than from a single 
pseudospherical surface (as considered in [5]). Denote 

m sin(w — Co) sin(cu + d;) 

n = r (1) - r = — —. — ^r f H —. — ^r„. (17) 

sin(2w) 5 sin(2w) v ; 

Then n is a unit vector tangent to both surfaces r and and determines what is 
called a pseudospherical congruence. Normal surfaces of this congruence are the constant 
astigmatism surfaces sought. 

Proposition 3. Let rj, c) be a general solution of system (14), where c is an 

integration constant. Then f = r — fh, where f = h^da/^/dc) and n is the unit vector 
given by formula (17), is a surface of constant astigmatism having surfaces r and 
as evolutes. 

Proof. The surface r = r — fh is normal to the congruence determined by the surface 
r and vectors n, if and only if f £ ■ n = f v ■ n = 0. By virtue of the Gauss- Weingarten 
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system (8) above, the derivatives of r = r — fn can be written as 

/. sin 2u (l) + sin2w „ sinfa/ 1 ) — oj) „ \ 
V 2 sin 2a; sin 2a; / 

/sin2(w (1) +w) , sin(w (1) +w 



sin 2a; / 



\ 2 sin 2a; ' sin 2a; 

sin — sin 2a; sin^ 1 ) + o>) , 
r,/ ~ M ? 2 sin 2a; _/ sin 2a; / "^ r " 



/sin^'-w sin^'-wj \ . (1) 

+ ( / /Jr^ + snuV J - u)fn, 

V 2 sin 2a; sin 2a; / 



>in2(o; ( ' — lu) sin(c 
2 sin 2a; ' sin 2a; 

Now it is straightforward to check that the conditions f ^ - n = and r v ■ n = reduce 
to equations (15). Moreover, it is a routine to verify that f is of constant astigmatism 
equal to 1. Actually an equivalent computation will be done in the proof of the next 
proposition under the same assumptions. □ 

Based on Bianchi's observation above, Proposition 3 shows that the constant 
astigmatism surfaces f = r — /n can be found by purely algebraic manipulations and 
differentiation once a one-parameter family of pseudopotentials is known. 

However, since £, 77 need not be curvature coordinates on the constant astigmatism 
surfaces f = r — /n, Proposition 3, as it stands, yields neither a solution of the constant 
astigmatism equation nor an orthogonal equiareal pattern on the sphere n. Yet the 
coordinates £,77 have a geometric meaning of a slip line field according to Definition 3. 

Proposition 4. Let u/ 1 '^, 77, c) be a general solution of system (14), where c is an 
integration constant, let f = k^du/^/dc) and x = df / dc. Let 77) be a solution of 
the system 

y% = e~f sin(o; + a/ 1 * 1 ), y v = sin(o; — a/ 1 - 1 ). (18) 

Then x,y are adapted curvature coordinates on the surface r. Moreover, if z = e~ 2 ? , 
then z(x,y) is a solution of the constant astigmatism equation (1). Finally, zdx 2 + dy 2 /z 
is an orthogonal equiareal pattern on the unit sphere h, while £,77 is the associated slip 
line field according to Definition 3. 

Proof. Continuing the routine computations started in the proof of Proposition 3 we 
obtain 

I = |(l + 2/ + 2/ 2 )(l - cos 2(w + u/ 1) ))d£ 2 
+ (1 + 2/)(cos2u;-cos2oj {1) )d£d77 
+ |(l + 2/ + 2/ 2 )(l - cos 2{u - a/ 1 ))) d£ 2 , 

fl = i(l + 2/)(l-cos2(a; + a; (1) ))d£ 2 
+ (cos 2lu — cos 2w (1) )d£d77 
+ Hl + 2/)(l -cos2(o;-ojW))d£ 2 . 
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The corresponding shape operator is 

/ , . l-cos2(w {1) -cj)\ 



1 

2/(7 + 1) 



2/ + 1 



cos 2a/ 1 ) — cos2w 



1- cos 2(^+0;) 2f + 1 

\ cos 2a;'- 1 - 1 — cos 2a; / 



Its eigenvalues are the curvatures, namely 1/ / and l/(/ + 1). We choose p\ — f + 1, 
P2 = / to have pi — p 2 = 1- The eigenvectors yield two principal directions 

n± = (cos2w (1) - cos2w)— ± (1 - cos2(w + w {1) )) — . 

It is easy to check that x and y are integrals of the fields 11+ and Il_, respectively. This 
implies that x, y are curvature coordinates. Now, equations 

I = u 2 dx 2 + v 2 dy 2 , II = (u 2 / Pl ) dx 2 + (v 2 /p 2 ) dy 2 . 

yield 



u 



(f + l)/e f , v = fe f . (19) 



Substituting into condition (2), we see that the curvature coordinates x, y are adapted. 
Moreover, it is straightforward to verify that z = e~ 2f satisfies the equation of constant 
astigmatism (1) with respect to independent variables x,y. 

To prove that £, rj is the associated slip line field on the sphere n, it suffices to check 

that 

cos3w — cosw — cos^u/ 1 * 1 — uj) + co&{2u^ + uj) 



nt = 



'< sin 2 uj cos uj 



cos(2w (1) + uj) - cos(2w (1) + 3uj) . , m s 
+ V — 4, V -r v + sin(w (1) + uj) n 



cos 3uj — cos uj + cos(2a; ( - 1 - ) — uj) — cos(2u;^ + uj) 
8 sin uj cos uj 

cos(2u; (1) - uj) - cos(2w (1) - 3w) , m 

H ^ " r £ - sin(w (1) -u)n 

8 sin cucosw 

bisect the right angle between 

cos(o; ( - 1 - ) — uj) cos(w < - 1 - ) + w) _ f 

n, = e — r f H e — r„ and n„ = e J n, 

sin2we / C sin 2uje f v y 

according to Definition 3. This is straightforward. □ 

Corollary 1. If S is a constant astigmatism surface, then the asymptotic coordinates 
on the focal surfaces of S correspond to slip line fields on the Gaussian image of S. 

Proposition 4 allows us to construct one of the adapted curvature coordinates by 
purely algebraic manipulations and differentiation, while the other curvature coordinate 
has to be obtained by integration. It is therefore natural to ask whether one could obtain 
superposition formulas for f,x,y similar to formula (12). The answer is positive. 
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Definition 4. Given two sine-Gordon solutions to and w^' related by the Backlund 
transformation B^ x \ let f( x \ x^ x \ denote functions satisfying the compatible 
equations 

/W = A cos(u^ + u) , 4 A) = i cos(w (A) - u) , 

4 A) = Ae' (A) sm(u^ + w ), = ^ (A) sm(oj^ - w ), (20) 



= Ae • / ' <A) sin(u/ A ) = -|e ^ (A) sin(w v " - x). 



(A) 



The quantities f( x \ x^ x \ will be called associated potentials corresponding to the 
pair u, uj( x \ 

Proposition 5. Let ui, u>^ Xl \ uj( X2 \ uj^ XiX ^ be four sine-Gordon solutions related by the 
Bianchi superposition principle (12). Then the associated potentials f( XlX2 \ x( XlX2 \ 
y(AiA 2 ) corresponding to the pair u^ Xl \ lo^ XiX ^ are related to the associated potentials 
f^ X2 \ x( x ' 2 \ y^ corresponding to the pair 00,00^ by formulas 

/ (A 1 A 2 ) = /( A 2 ) _ ln ^cosM Al ) -U^) -^-^), 

(AlA2) _ A:A 2 / {A2) _ 2A 1 A 2 sinM A i>-c^) /( , 2 A 



A? - Xl V A? - 2AiA 2 cos(co^ - u (x ^) + \\ 

y (A!A 2 ) = (^1 _ ^U^) - 2e^ (A2) sin(a/ Al > - W ( A2 )), 
VA 2 \\' 

up to an additive constant. 

Proof. It is straightforward to check that f( XlX2 \ x^ XlX2 \ y( XlX2 *> given by formulas (21) 
satisfy 

/( AlA2 ) = A 2 cos(J x ^ + uj^), 6 X ^ = ±- cos(a/ AlAa > - u/ Al >), 

A 2 

x pW = A 2 e' (A) sm(u^ + u ™), x^ = l e ^ (AlA2) sm(J x ^ - u^), 

A 2 

y ^ = A 2 e^ (AlA2) sin(u/ AlA2 ) + u^), y^ = -le"^ 2 ' s in(J x ^ - J x ^). 

A 2 

whenever f( X2 \ x^ X2 \ y^ satisfy (20) with A = A 2 . □ 

Example 3. (One-soliton solutions) Let us apply the procedure outlined above to the 
one-soliton solutions 

W (A) = #<A)(o) = 2arctanexpp A , 

of the sine-Gordon equation. Here and in what follows we denote 

Pa = A£ + -, q = £-r}. 
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As is well known, these one-soliton solutions correspond to the Dini surfaces (helicoids 
of the tractrix) 

2A /sechp A sing\ ^ _ ^ 2 /0 N 
r (A) = - - ^ 2 sechp A cosg + 
\px - tanhp A / 

We now proceed to the complementary surfaces of the Dini surfaces, which correspond 
to the nonlinear superposition of ui^ and wW, i.e., the two-soliton solutions 



w (Ai) = J B(i)( w (A)) = 2arctan 



(A + l)(e PA 



c 



(A-l)(l + e^ +Pl+c )' 

where, obviously, pi = £ + ?y. The particular case of A = 1 (the Beltrami pseudosphere) 
is excluded from consideration. 

After tedious computations, one obtains the resulting quantities x = x^ xi \ y = y^ xi \ 
z = e~ 2 f 1 . They are 

A (A - \) 2 {c 2 A 2 B 2 - cQ - (A + lf{ Cl B 2 + c 2 A 2 ) + 4( Cl - c 2 )AAB 
X ~A 2 -1 X (\-l) 2 (A 2 B 2 + l) + (\ + l) 2 (B 2 + A 2 )-8XAB ' 

_ 41n5 (A 2 + l)lnA AX(AB + 1)(A - B) + c 3 (ci + c 2 )(A 2 - 1)4(1 + B 2 ) 
V ~ c 1 + c 2 ( Cl + c 2 )A + ( Cl + c 2 )Ae s (l + e 2 ^) 

_ / (A - 1) 2 (A 2 E 2 + l) + (A + 1) 2 (£ 2 + A 2 ) - 8AA£? \ 2 
* V ( Cl + c 2 )AA(l + J B 2 ) " J' 

where A = e Pl = e^ +v and S = e PA = e x ^ +v ^ x , while ci, c 2 , c 3 are arbitrary constants. By 
eliminating £, 77 one obtains 

1 /4(AB + 1)(A-S) n A 2 + l. , _\ (A 2 -l)c 3 



y = — 

Cl 



where 



A(A 2 + l)( Cl + c 2 )y^ - y/fc 
(A 2 -l) 2 + (A 2 :r-Ac 2 -:r)V 

B = 2A 2 (c 1 + c 2 )y / ^ + y/fc 

~ (A 2 - l) 2 + (A 2 x - Ac 2 - x)(\ 2 x + Aci - x)z 

k = -[(A 2 - l) 2 + 2(ci + c 2 )A 2 v^ + (A 2 x - Ac 2 - x)(A 2 a; + Ac x - x)z] 

x [(A 2 - l) 2 - 2(ci + c 2 )A 2 \/^ + (A 2 a; - Ac 2 - x)(A 2 x + Aci - a;)*]. 

Eq. (22) is an implicit formula for a solution z(x, y) of the constant astigmatism equation. 
Using Proposition 3 it is now easy to construct the surface of constant astigmatism from 
its two evolutes r^ A ^ and r^ A1 ^ as well as the orthogonal slip line net on the Gaussian 
sphere n = r^ A1 ^ — r^, part of which can be seen on Fig. 4 (actually, the sphere 
is multiply covered). This example also demonstrates that sphere's slip line fields are 
prone to developing singularities. 

Construction of constant astigmatism surfaces related to the sine-Gordon n-soliton 
solutions is postponed to a separate paper. 
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Figure 4. Sphere's slip line field with features 



5. Lipschitz surfaces in principal coordinates 

In 1887 Lipschitz [22] presented a class of surfaces of constant astigmatism in terms of 
spherical coordinates related to the Gaussian image. To find the corresponding solutions 
of the constant astigmatism equation, one has to obtain the surfaces in terms of the 
pricipal coordinates. Redoing the computation is easier than transforming the Lipschitz 
result. 

Consider the unit sphere n = (cos sin 9, sin sin 9, cos 9) parameterised by the 
latitude 9 and longitude <fi. To specify an orthogonal equiareal pattern we let 9, <fi denote 
yet unknown functions of parameters x, y. Lipschitz defines a Stellungswinkel to be the 
angle u between n e and n x = <p x n^, + 9 x n e . The Lipschitz class is specified by allowing 
the Stellungswinkel to depend solely on the latitude 9. 

Theorem 1. The general Lipschitz solution of the constant astigmatism equation (1) 
depends on four constants hn, hio, h i, /i o and consists of functions 

where 

h = h n xy + h w x + h 0l y + h 00 , Hi = h u , H 2 = h u h 00 - hi h 0l . 

The constants H\,H2 are invariant with respect to the translations % x ,1. y , the scaling 
& and the involution 3. Formula (23) covers all Lipshitz solutions except a particular 
solution 

1 

Ci - \x - C J 
Ci,c being arbitrary constants. 
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Proof. Computing the first fundamental form in two ways we have 
d9 2 + sin 2 9 d0 2 = z dx 2 + (1/z) dy 2 , 

i.e., 

e 2 x + 2 sin 2 e = z, e x e y + (f) x (f) y sm 2 e = o, e 2 y + 2 sin 2 e = i/z. 

Eliminating z, we have 

[9 2 x + 2 sin 2 6){6 2 y + 2 sin 2 9) = 1, 9 x 9 y + 0,0, sin 2 9 = 0. (24) 

Since 

= (cos cos 9, sin cos 0, — sin#), = (— sin sin 0, cos sin 6>, 0), 
ria; = (9 X cos cos 9 — <p x sin sin X sin cos + 0^ cos sin 6>, — X sin#), 
we easily compute 

■ 9 X 

cosw = - — n — r = — i (25) 

Honce, the Lipschitz' condition amounts to z /6 l :r being a function of # alone. This can 
be conveniently written in the form 

0, = —J x (26) 
sm9 

(intentionally leaving sin# unabsorbed). Conditions (24) and (26) combine into the 
system 

0x= shi/ x ' <Py = ~(sm9)(Q 2 + lW ° y= (Q 2 ll)9 x (2?) 

Computing the compatibility conditions for we get (sin 9)9 XX + (cos#)# 2 = with the 
general solution 

9 = arccos h, h = h\X + h , 

hi, h being arbitrary functions of y. Now the third equation of (27) implies 

6 2 + 1 y 1 + (h lX + h ) 2 \dy dy J K ' 

Since the left-hand side is a function of 9 alone, the same is true for the right-hand side. 
Checking the Jacobian, we get 

hi /dh\ dh 2 \ 

0_ (i+(h lX +ho) 2 r/ 2 w x+ w)' { ] 

If hi = 0, then h = h is a function of y alone and from system (27) we easily obtain 
ho = h 01 y + h o and z = l/^i — (y + h 00 /hoi) 2 , which is the particular solution. 
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Otherwise hi ^ 0; it follows from (29) that hi, ho are linear in y and we can write 

h = huxy + h w x + h iy + h 00 , (30) 

with hij being arbitrary constants. Now, eq. (28) gives 

G> = Hih H 2 ^ where Hi = h H ^ = / ln / l00 _ h w h m . 
B + 1 1 — h 

Hence, 

1 - h 2 + V(l-h 2 ) 2 -4(Hih-H 2 ) 2 



B 



2(Hih - H 2 ) 



Since Hi, H 2 are constants, B is a function of h and, consequently, of 6 alone. Inserting 
into system (27) and then into z = 9 2 + 4> 2 X sin 2 9, we obtain the general solution (23). 

□ 

It is now easy to obtain the corresponding orthogonal equiareal pattern. 

Theorem 2. The orthogonal equiareal pattern corresponding to the general Lipschitz 
solution is n = (cos <fi sin 9, sin <fi sin 9, cos 6), where 

9 = arccos h, 

\n(h n x + hoi) , f 1 - h 2 + V(l - h 2 ) 2 - 4(Hih - H 2 f , . 

= Vn + J 2(Hih-H 2 )(l-h 2 ) d '' (31) 

h = hnxy + hi x + h iy + h 00 , H x = hu, H 2 = hnh 00 - h 10 h 01 . 

Proof. We need to know 0, i.e., we have to integrate the first two equations (27). It is 
easily observed that 4> x h y — 4> y h x + 1 = 0. Solving this PDE for 0, we obtain 

4> = - Hhll * + hoi) +*(h), 

h u 

while for $(/i) we get 

d$ _ 1 - h 2 + y/(l - h 2 ) 2 - A(Hih - H 2 ) 2 
dh ~ 2(Hih - H 2 )(l -h 2 ) 

□ 

The Stellungswinkel u is a function of the lattitude 9 as required; namely 
1 1 - h 2 + y/(l - h 2 ) 2 - A(Hih - H 2 ) 2 



cos 2 u 



B 2 + 1 " 2(1 - h 2 ) 

sin 2 9 + V / sin 4 ^-4(iJ 1 cos^- J ?J 2 ) 2 



2sin 2 ^ 
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Remark 2. It is easy to check that the general Lipschitz solution (23) satisfies 

h u s + h 01 i x - h 10 t y = 0, 

where t x = z x , t y = z y , s = xz x — yz y + 2z are generators (see, e.g., [10]) of the 
Lie symmetries % x , 1. y , &, respectively. This means that (23) is a symmetry-invariant 
solution of the constant astigmatism equation. 

Example 4. When the integral in (31) can be expressed in terms of elementary 
functions? Assuming that hn is nonzero, hio and hoi can be removed by shifts, so we set 
hio — hoi = 0. Consider the expression under the square root in (31). Its discriminant 
with respect to h is proportional to 

(1 + Hi + 2H 2 ){\ + Hi - 2# 2 )(# 1 - # 2 ) 2 (#i + H 2 f 

= h 4 n (l + h 2 n + 2/m/ioo)(l + h 2 n - 2/in/ioo)(l - h 00 ) 2 (l + /ioo) 2 , 
which is zero if and only if 

1 + h 2 

h 00 = ±1 or h 00 = ±- 



2h 



ii 



In these cases, <fi can be expressed in terms of elementary functions. For hoo = ±1 we 
have 

VT^C 2 , A(l - 2C 2 ± h + v / T 3 C^ T I) 2 - 4(C 2 T h) ) ' 
— -F m 



2C \ hTl 



ln(Cx) ln[/i 2 - 1 + (h T 1) V{h T l) 2 - 4(C 2 =f /i) ] 
C + " ~^2C 



1 fV{h±l) 2 -AC 2 
±-arctan^ 

where C — H ± — h u is a constant and h = Cxy + 1. The orthogonal equiareal pattern 
corresponding to hoo = 1 and /in = 1/4 can be seen on Fig. 5. 
In the second case, when h 00 = ±(1 + hl 1 )/2h n , we obtain 

ln(Cx) \n(h±C+ \/h 2 ±2Ch-C 2 -2) 
r C 2C 

1 , f-C 3 - 3C ± /i(3C 2 + 1) - (1 - C 2 ) \A 2 ± 2Cfc - C 2 - 2 
In 



2C V (C 2 =f 2C7i + 1) 



1 ( 2(Ch ± 1) y/h 2 ± 2C7i -C 2 -2 

± 2 arCtan V /> 2 - - C 2 ± 4C7* - 3 

Here h, C have the same meaning as above. No figure is provided in this case, since 
and 6 cannot be simultaneously real. 
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Figure 5. The orthogonal equiareal pattern on the sphere corresponding to one of the 
Lipschitz solutions 

6. Conclusions 

Summarizing, we identified the constant astigmatism equation and the sine-Gordon 
equation as integrable models in two-dimensional plasticity on the sphere with respect 
to parameterisation by principal stress lines and slip lines, respectively. We remark in 
this context that the majority of exact solutions in plasticity that can be found in the 
literature are either due to linearisable systems or come from symmetry methods; see, 
e.g., [3, 20] and references therein. 

We have also extended the classical Bianchi superposition principle so as to be able 
to generate solutions of the constant astigmatism equation by algebraic manipulations. 
Finally, revisiting the classical Lipschitz surfaces of constant astigmatism, we have 
identified them as corresponding to invariant solutions. 

In conclusion, we are able to say that obtaining large families of exact solutions 
of the constant astigmatism equation as well as interpreting them as plastic flows on 
a sphere is merely a matter of routine. Since computations quickly leave the realm of 
elementary functions, the examples and illustrations scattered throughout this paper 
are only the simplest ones. More are to follow in a subsequent paper. 
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